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Sets of orthogonal basis functions over circular areas — often representing pupils in optical 
applications — are known in the literature for the full circle (Zernike or Jacobi polynomials) and 
the annulus. Here, an orthogonal set is established if the area is two non-overlapping circles of equal 
size. The main free geometric parameter is the ratio of the pupil radii over the distance between 
both circles. Increasingly higher order aberrations — as defined for a virtual larger pupil in which 
both pupils are embedded — are fed into a Gram-Schmidt orthogonalization to distill one unique set 
of basis functions. The key effort is to work out the overlap integrals between a full set of primi- 
tive basis functions, which are chosen to be products of powers of the distance from the mid-point 
between both pupils by azimuthal functions of the Fourier type. 
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I. AIM AND SCOPE 

Manufacturing schemes of lenses and mirrors inevitably prefer circular cross-sections of beams, and the associated 
description of functions (aberrations) defined across these fields calls for basis functions on this circular support, the 
best- noted probably being the Zernike functions [1, 0, [13, [H, HH . Masking a central circular portion of a circular 
beam leads to annular, ring-shaped regions, for which orthogonal basis sets are also established fsl. fiol. [l^. [2^ [26| . 

This work proceeds to the task of defining such a basis set for a two-beam interferometer, in which the input pupil 
is defined by two disconnected circular areas of equal radius [ll| . 

We define the area of integration in a global spherical coordinate system centered in between the two apertures 
in Section [TTl Supposed anonymous functions defined over these apertures are expanded with a separation ansatz as 
products of powers of the distance to the origin of coordinates by the usual Fourier scries in the azimuth, all integrals 
over products of these can be reduced to a generic integral, summarized in Section Hill The value of this article lies in 
the the reduction formulas of two associated integrals in two appendices. Section HVl proceeds with an application, the 
re-orthogonalization of the Zernike basis functions — defined over the larger area that encompasses both circles — with 
respect to the two circular regions that define the binocular pupil. 

II. BINOCULAR GEOMETRY 

Orthogonality of functions fk and // over two-dimensional areas is defined through their product integrated over 
the area 

{fkJi) = J J rkfidn^Ski- (1) 

In Cartesian coordinates x and y, or circular coordinates with distance p to the origin and azimuth 9, 

X = pcosO; y — psinO, (2) 

the differential is dil = dxdy or d^l = pdpdO. This manuscript deals with two-dimensional areas that are the sum of 
the interior of two circular pupils represented by 

\2 I „.2 ^ 2 



{x±RY + y' <t', (3) 
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FIG. 1: The two-dimensional area and the Cartesian and polar coordinates covered by the two circular sub-pupils. The center 
of symmetry and origin of the global circular coordinates is at x = ?/ = 0. In the first quadrant, the limit 6 < arcsing of 
integration is given by the tangent of a radial vector on the pupil with index P = -1-1 along the dotted line. For general 9, the 
dashed line illustrates how the radial vector intersects the front and rear side of the circle, associated with the two signs in (f6]l. 



where R is half the distance between the two pupil centers, where 2R is the interferometric baseline, and r is each 
pupil's radius (Fig. [T]). This transforms the integral operator into a sum over both circles. 



R+r 



pdp 



R-r 



de- 



2pflcose<r2-_R2_p2 



de 



2p_Rcose>p2-r2+i?2 



(4) 



Scaling distances in units of i?, z = p/R, leaves one essential shape parameter, q = r/R. To avoid double-counting of 
areas, these must not overlap: 

0<q<l. (5) 
The range of radial distances that lie inside the circle centered at x = +R for some fixed direction 9 is 



z — cos 9 ± \ — sin 9, 



and for the other one centered sX x = —R, 



z = — cos 9 ±\J q^ — sin 



(6) 



(7) 



In each of the two equations, the lower sign connects to the intersection with the circle rim that is closer to the origin, 
the upper sign to the intersection with the farther one. 



III. GENERIC INTEGRAL 



At the heart of this work is performing overlap integrals ([T]) in analytical terms over the area described above; since 
this will be based on spanning the functional space with "primitive" basis functions of the type 

Pfc=p„,,„,(z,0) = z"'=e™'=^ (8) 

the generic integral reads 

Ar,.ra{q)= [ d9f zdzz^e'^\ n, m = 0, 1, 2, . . . (9) 

/i . n\ ^ I2zcos0<q —z —1 

J sill 6 \<.q J ^ Q Q 

' 2zcose>z^-q^+l 
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i is the imaginary unit, m the azimuthal frequency, and n the power to build a complete radial basis. With a sign 
tag P defined for both circles. 



P = 



+1, \0\<7r/2, 
-1, 7r/2 < \e\ < IT, 



(10) 



the integration over 2"+ may be executed, 



(n + 2)An,,niq) = / de 

'\sine\<q 



Pcos0 



q"^ — sm 



n+2 

- (Pcos6l - 



^2 - sin^ ( 



n+2 



(11) 



The contribution from the imaginary part proportional to s'm{m6) vanishes; e*™^ can be replaced by cos{m9) in this 
equation. Considering the coordinate transformation y <-> —y or —9, the integral over all four quadrants can be 
reduced to an integral over the first and second quadrant and a factor of 2: 



{n + 2)An,miq) 



2 1 de 

sin ^? I < g 
0<9<7r 



Pcos9 



n+2 

g2-sin^6l) - (Fcos( 



— sin^ I 



n+2- 



cos{m9). (12) 



Considering also the variable transformation 9 ^ n — 9, the parities of the Chebyshev term Tm(cos6') = cos{m9), of 
cos 9 and of P, this vanishes for odd m and reduces to an integral over the first quadrant for even m. 



(n + 2)A„,„((?) 



sin 9<q 
0<e<ir/2 



d9 



cos ( 



— sin^ I 



n+2 



— I COS W — \/ — sin^ i 



n+2- 



cos{m9). (13) 



If we define the integrals 



B 



n+2,m 



(q) 



J d9\cos9±\Jq^ ~ sin^ 9j T„, (cos 



(14) 



this can be rephrased 



(n + 2)A„,„,(g) _ f ,m= 1,3,5,7,... 

4 \ K+2,m{q) - K+2Aq) ,m = 0,2,4,6,... ■ 

Table [J summarizes the results obtained in Appendix VK\ as series expansions of ^„,m('?) for small n and m. 



(15) 



n m An,,n(g)/7r 
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+2q" 
+2g2 - 

+2g2 - 4g* + 2g« 

+2g2 - Qq^ + 12g® - 5g* 

+2g^ - IGg" + 40g*^ - 40g* + 14g^'' 

+2g2 - 25g'' + lOOg'' - ll^q^ + 140g"' - 42g^2 



49/196608g^ 



+2g^ + l/4g'' + l/96g'' + l/512g'* + 5/8192g"' 
+63/524288g" + 1089/16777216g^'^ + . . . 
+2g2 - 3/4g'' + 3/32g® + 5/512g'* + 21/8192g"' + 63/65536g^^ 
+231/524288g" + 3861/16777216g^'' + . . . 

+2g2 - 15/4g'' + 75/32g^ - 175/512g* - 315/8192g"' - 693/65536g^2 

-2145/524288g" - 32175/16777216g^® + . . . 

+2g^ - 35/4g^ + 1225/96g® - 3675/512g'* + 8085/8192g^° 

+7007/65536g^^ + 15015/524288g" + 182325/16777216g^^ + . . . 

+2g2 + g4 

+2g2 

+2g2 - 3g'' + 2g'' - l/2g* 
+2g2-8g^ + 12 
+2g2 - ISg'' + 40g'* - 50g* 

+2g2 - 24g'' + lOOg'' - 200g** + 210g^'' - 112g^^ + 24g^* 
+2g2 + 9/4g'' + 3/32g'' + l/512g'* + 9/40960g^" + 3/65536g^^ 
+7/524288g" + 81/16777216g" + . . . 

+2g2 + 5/4g'' - 5/32g'' + 5/512g'* + 5/8192g"' + 7/65536g^^ 



8g» + 2gio 



30gi° - 7gi2 
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+15/524288(7" + 165/16777216?^'^ + . . . 
3 4 +2q2 - 7/4?'' + 35/32g'' - 175/512q* + 245/8192g^" + 147/65536g^^ 
+231/524288g" + 2145/16777216g^'* + . . . 

3 6 +2g2 - 27/49"* + 315/329** - 3675/512g* + 19845/8192gi° 

-14553/65536?^^ - 9009/524288?" - 57915/16777216?^® + . . . 

4 +2?2 + 4?* + 2/3?® 
4 2 +2g2 + 3?"* 

4 4 +2g2 

4 6 +2?2-5?'* + 20/3?''-5g* + 2gi°-l/3g*^ 

4 8 +2g2 - 12?* + 30?® - 40?* + 30?^° - 12?^^ + 2?" 

4 10 +2g2 - 21g* + 84g® - 175g* + 210?*" - 147g^^ + 56?" - 9g^® 

5 +2g2 + 25/4g'* + 75/32?*' + 25/512g** + 5/8192g"' + 3/65536g*^ 

+25/3670016g" + 25/16777216g*'' + . . . 
5 2 +2?2 +21/4?* + 35/32?®- 35/512?® +21/8192?*° + 7/65536?^^ 

+7/524288?" + 45/16777216?*® + . . . 
5 4 +2g2 + 9/4g* - 21/32g® + 105/512g® - 315/8192?*° + 147/65536?*^ 

+63/524288?" + 297/16777216g*® + . . . 

5 6 +2g2 - ll/4g'* + 99/32?® - 1155/512?® + 8085/8192?*° 

-14553/65536?*^ + 7623/524288?" + 14157/16777216?*® + . . . 

6 +2?2 + 9?* + 6?® + 1/2?® 
6 2 +2?2 + 8?* + 4?® 

6 4 +2?2 + 5?* 
6 6 +2?2 

6 8 +2?2-7?* + 14?®-35/2?® + 14?*°-7?*2+2?"-l/4?*® 

6 10 +2?^ - 16?* + 56?® - 112?® + 140?*° - 112?*^ + 56?** - 16?*® 

7 +2?2 + 49/4?* + 1225/96?® + 1225/512?® + 245/8192?*° 

+49/196608?*^ + 7/524288?" + 25/16777216?*® + . . . 
7 2 +2?2 + 45/4?* + 315/32?® + 525/512?® - 315/8192?*° + 63/65536?*^ 

+15/524288?** + 45/16777216?*® + . . . 
7 4 +2?2 + 33/4?* +99/32?® -231/512?® +693/8192?*° -693/65536?*^ 

+231/524288?** + 297/16777216?*® + . . . 

7 6 +2?2 + 13/4g* - 143/96?® + 429/512?® - 3003/8192?*° 

+7007/65536?*^ - 9009/524288?** + 14157/16777216?*® + . . . 

8 +2?^ + 16?* + 24?® + 8?® + 2/5?*° 
8 2 +2?2 + 15?* + 20?® + 5?® 

8 4 +2?2 + 12?* + 10?® 
8 6 +2?2 + 7?* 
8 8 +2?^ 

8 10 +2?2 - 9?* + 24?® - 42?® + 252/5?*° - 42?*^ + 24?** - 9?*® 

9 +2?^ + 81/4?* + 1323/32?® + 11025/512?® + 19845/8192?*° 

+1323/65536?*^ + 63/524288?" + 81/16777216?*® + . . . 
9 2 +2?2 + 77/4?* + 1155/32?® + 8085/512?® + 8085/8192?*° 

-1617/65536?*^ + 231/524288?" + 165/16777216?*® + . . . 
9 4 +2?2 + 65/4?* + 715/32?® + 2145/512?® - 3003/8192?*° 

+3003/65536?*^ - 2145/524288?" + 2145/16777216?*® + . . . 

9 6 +2?2 + 45/4?* + 195/32?® - 715/512?® + 3861/8192?*° 

-9009/65536?*^ + 15015/524288?** - 57915/16777216?*® + . . . 

10 +2?^ + 25?* + 200/3?® + 50?® + 10?*° + 1/3?*^ 
10 2 +2g2 + 24?* + 60?® + 40?® + 6?*° 

10 4 +2?^ + 21?* + 42?® +35/2?® 
10 6 +2?2 + 16?* + 56/3?® 
10 8 +2g2 + 9g* 
10 10 +2?^ 

TABLE I: Series expansion of jAn^rn (?)/7r. It is a polynomial in ?^ if 
n + m is even, else shown in truncated form up to 0(?*®). 



To lowest order in Afi^jYi equals 27r?'^, the total area of both circles. This is expected, because for circles far 
away from their common origin at ^ = 0, the relative strength of the variation introduced by the power n and by the 
modulation oc cos{m0) of the function in the integral kernel loses importance. 
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IV. GRAM-SCHMIDT ORTHOGONALIZATION 



A. Procedure. 



The integral evaluated in Section IIIII allows to calculate the overlap integral (inner product) between any two 
functions expressed as linear combinations ( "contractions" ) of "primitive" basis functions of the form ^ in the global 
circular coordinate system centered at the middle between the two sub-pupils, because the overlap between two of 
these is 



The quickest, obvious way of obtaining some orthogonal basis set from any set of contracted primitive basis function 
is to diagonalize the overlap matrix containing all the overlap integrals between pairs of these basis functions 0- To 
end up with some standardization of these orthogonal bases, we use the Gram-Schmidt procedure, which builds this 
set {/fc} incrementally. At each step of the procedure, an ansatz 



fk+iiz,d) = Pk+l 



gk+i{z,e) + ^jk+i,ifi{z,9) 



1=1 



(17) 



is made for the next, (A: -I- l)st additional basis function fk+i, given a seed function gk+i{z, 9) plus the /; generated by 
the earlier steps. Essentially, the projections of the seed along all earlier directions are subtracted, and the residual 
is normalized to unity. The request of orthogonality 



f*fidn = {fkji) = 5ku k,i = i,...k + i, 



(18) 



means the projection coefficients 7 can be calculated from the k overlaps between the seed and the earlier basis 
functions, 

lk+i,i = "{gk+i.Ji), l = l,...k. (19) 
The normalization /3fc+i is finally computed from the self-overlap of the seed and the sum over the squared 7, 

k 



1 = Pl+1 



{9k+i,gk^ 



"^^k+i.. 



(20) 



The Gram-Schmidt methodology is well known in the literature for different geometric shapes of pupils [15|, |23|, |25| . 

Still, the procedure establishes different basis sets depending on the order in which the seeds g are fed into the 
procedure, and depending on which functional subspaces they span. 



B. Zernike Seeds. 



To define a unique set {fj} of functions orthogonal over the binocular pupil, we may choose the real and imaginary 
parts of the primitives ([8]) in increasing order of complexity, ie, increasing order of aberration and increasing n and to, 
as the seeds, as if one would subduce the Zernike polynomials Zj(z, 9) over the full super-pupil of radius R (normalized 
to z = 1) in the Noll order of indexing [l3| into the two sub-apertures. (We identify the variable z with the radial 
variable of the Zernike polynomials, although the corresponding Zernike radius would need to be 1 -I- g, not 1, to 
cover both sub-apertures in full.) There is an "outer" loop over n = 0,1,2,... and an "inner" loop over m = n 
(mod 2), . . . n, considering only even n — m: 

51 = 1; g2=zcos9; 53 = z sin 6*; (21) 
34 = 22^-1; = z2sin(20); ge = z^ 008(29); (22) 
gT = {3z^ -2z)siii9; gs = {3z^ ~ 2z) cos9; gg = z^ sm{'S9): gio = z^ cos{W); .... (23) 

For this choice, the n + 1 basis functions from gi+„(„+i)/2 up to and including g(n+i){n+2)/2 ^-re associated with a 
polynomial of order n in z. The arithmetic remains real- valued, because the "atoms" of the seeds are the separated 
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real and imaginary part of ((T6)) is split into 



z"" sin(TOfc6')z"' sm{mie)dn ^ 
z"" cos(mfc6')z"' sm{mie)dn = 



-A 



nfc+n,,|?Tifc-m,| J 



(24) 
(25) 
(26) 



The first basis functions created with this recipe are discussed shortly in analytical form. We start with the "global 
common piston" gi which just needs to be normalized: 



/i 



1 1 



q V27r 



(27) 



Next we feed what represents most of the differential piston, g2, which turns out to be already orthogonal to /i and 
only needs to be normalized, 



-z cos 9. 



qV tt{4 + g2) 

Next we feed 173, some common sideways tilt perpendicular to the baseline between the two pupils, 

/s = \ \ -zsind. 
q'' \ TT 



(28) 



(29) 



The first case of nonzero overlap with an earlier basis function occurs when we use 34 as a seed (some nodding tilt 
between the sub-pupils along the baseline), which has a nonzero component along /i: 



1 



6 



g2V ^(12 + g2) 



z2- 



2 + g2 



Feeding 55 to gs we get 



q'^ y 7r(6 + q"^) 



z^ sin(26l). 



(30) 



(31) 



g V 7r(18 + g2)(i2 + q2) 



^^±^z'cos{2e)^^z' + 5 



(32) 



q-^ y 7r(12 + q'^) 



(33) 



q^y TTQs{4 + q^: 



[3(4 + q'^)z^ cose - {12 + 21q^ + 2q^)zcos9] . 



The abbreviation Qs = 288 



.4 I „6 



q is used. Continuing with 179, 



g2 y 7r(24 + g2)(-i2 + g2) 



12 + g2 



36 



z-^ sm{30) + 21zsin(6l) - -rrz sin 



(34) 



(35) 



'10 



q V 7rQ8(2304 + 704q2 + 480g4 + 60q6 + q») 



^z^ cos{3e) + 21{4: + 2q + q^){4:-2q + q^)zcos{e) - 



12(24 + 8(72 + 3g4 



-z cos f 



(36) 
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/ii = 



10 



ttQh (18 + g2) 

X [3(18 + q^)z'^ - 3(24 + 2&q^ + q')z^ + 
The abbreviation Qn = 2880 + 240?^ + 78g4 + g6 jg ^^^^ 



4 2 , 108 + 30g2 + 369^^ + 



- 12(3 + g2)z2cos(26i) . (37) 



/i 



9^ V ttQh (64000 + 10240g2 + 1840^4 + 128q6 + gS) 



4Q 



11 _4 



cos(26') 



3(3840 + 3200g2 + 4005^^ + 84g6 + q») ^ 



cos(26i) - 21(160 + 108g^ + 4g* + q") 
r 

12(960 + 1040g2 + 38^^ + 13g^) _2 60(192 + 8^^ + 3^^) _4 



(38) 
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q^ Y 7r(640 + 80g2 + 56g4 + q^){6 + q^) 
Figs. HHH illustrate the first eighteen of these functions for q = 1/2 



4(6 + q^)z^ sin{2e) - 3(8 + 12^2 + q'^)z'^ sui{2e) . (39) 



INTERFEROMETRIC SIGNAL 



Expansions in orthogonal bases lead to accelerated book-keeping: the integral over the square of a function becomes 
the sum of the squared expansion coefficients (Parseval's equation). In Maxwellian electrodynamics, the function is 
one of the two polarizations of the electric field vector f , and the simplification addresses how much total energy passes 
through the cross-section. If the circles are the entrance to a two-beam interferometer, this addresses computation of 
the photometric signal. 

For a field expanded in the primitive basis, 



£iz,e)^J2 



(40) 



the interferometric signal correlates values at conjugated points P ~ ±1, x = P + s cos Lp, y = s sin 1^9, sharing the same 
local radial coordinate s and azimuth Lp (Fig.[T]). The coordinate transformation to global {z^6) circular coordinates 
are 



Vl + 2Pscos</7 + s2; 



arctan(y/a;) — (7r+) arctan 



s sin(/3 



P 



cos if 



(41) 



where (-I-tt) indicates that tt is to be added for P = —1 if arctan denotes the principal value. The interferometric 
signal (spatial autocorrelation) is calculated by multiplying £ at two conjugated points in the pupils, a distance 2R 
apart in the a;-direction, and integrating over s and (p. After transformation of £ at (z, 9) to the individual (s, ip) 
coordinates, the interferometric signal breaks down into a sum over products of the expansion coefficients in terms of 
these shifted/scaled polynomials. As we have set up each fk as a linear combination of Zernike Polynomials Zj<k{r, 0), 
the route to transformations to polar coordinates originating at the two circle's centers is known from the literature 
01 Hi [Hj HOj l22|]. We summarize this in our notation; this is off-topic in the sense that it is not related to the 
orthogonality introduced above. 

The interferometric intensity is J £p^^iis, ip)£p=+i{s, ip)sdsdip — ; rj'^rjilk^i, the generic information contained 



sds 
Jo 



d(/;(z"''e-""'=^)|P^_i(z"'e™'^)| 



(42) 



''\ds [ dip{l - 2sC0Sy) H- S^)"'°/^(l + 2sC0Sy) -f s2^r.,/2^-^rn,[^+^vct^n _ y, ] g^"^ ' ^^'^^^^ _ (43) 

3 "'0 
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FIG. 2: The basis functions /i to /e for q = 1/2. The level of /; = is indicated with small horse-shoe rims around both 
sub-pupils. 

We expand the integrand in a power series of s, integrate term by term and show the results in form of the first terms 
of a power series in q. The mean and excess of the four parameters, 

m(+) = (mfc + mi)/2- m^") = (wfe - m;)/2; = K + n;)/2; n^") = (n^ - n;)/2 (44) 




FIG. 3: The basis functions fr to /12 for q = 1/2. 



FIG. 4: The basis functions /13 to /is for q = 1/2. 
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are defined to compress the notation. 



nq 



2m(+) - n 



t V 



t 



,6 

192 



t nq' 



2m(+) 



,(-) 



8 + (2m(+) - + 6(2m(") - n^) 



9216 



+12(2m(") - n(+))(4 + 2m("' - n(+)) 



t 



nq 



10 



737280 



2m(+) 



,(-) 



384 + (2m(+' - n^-^f + 20(2™^ _ n("))2(4 + 2m^-^ - n^) 



88473600 



+20(2m(") - |20 + 3(2m(") - 
(2m(+) - n'^-^Y + 10(2m(+) - n'--yf |l6 + 3(2™^-) - nW)| 
+4(2m(+) - n^-^ f |736 + 420(2m(") - + 45(2m(~) - n^+^ f^ 



+120(2m(") - n(+))(8 + 2m'--^ - nW)(4 + 2m(") - 



t 



7rg 



14863564800 



(45) 



The symbol f indicates that bracket to its left represents the product of two factors. The first factor is the bracketed 
term as written; the second factor is the term after the substitutions n^") — > — n^^-* and m*-") — > — m*-^-*. So in the 
second factor of the product, some sign flips occur whenever the total power of the m^"^ and n'--~^ is odd. 

It is equivalent to two re-expansions of the field at shifted centers of the sub-pupils 6] followed by areal integration. 
Selection rules are implicit; if 2m*^^) — or 2to(+-' — n'--^^ vanish, many coefficients in the power series become zero. 



VI. FOURIER REPRESENTATION 



All basis functions fj are linear superpositions of the primitive type ([5]) after the cos(m0) are replaced by [exp(im6) + 
exp(— im0)]/2 and the sin(m0) by [exp(im0) — cxp(— zr7i0)]/(2i). The two-dimensional Fourier Transform of these is 



An.miq,cr) = J exp{2nicr ■ z)z"' exp{im9)dn 



(46) 



for some wave number cr. These integrals are calculated by individually translating each of the two circular areas to 
the origin of coordinates as described in (|¥l]) . 



P=±l s<q 



(47) 



where ax is the component of a along the baseline axis. The coordinate transformation may interpret ze as a 
complex variable. Since only the cases of even n — m need to be considered. 



is expanded in multinomials of x and y, x is replaced by a; -f P. Finally the substitutions 



y — s sm 



-tip] 



(48) 



(49) 



express z'^e^™ 



in the circular coordinates centered at x = P with radial coordinate s and azimuth ip. Even powers of 
P are dropped because P'^ = 1. Table HIl demonstrates the cases for small n and small non- negative m. 
This reduces each An^m{q, cr) to a flnite sum of integrals over centered circles of radius q, 



s<q Jo Jq 



(50) 



2^e™^'*l'"l J^' s"+ i2nsa)ds = 2^e""^'^l"l (2^)^+2 g»+i,l^l (2^^^), (51) 
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n m z" exp{imO) 

1 

1 1 P + se''^ 

2 2 Ps cos (ifi) + 1 + s'^ 

2 2 l + 2Pse''^ + s^e^"^ 

3 1 P + 2Ps'^ + Ps'^e'^''^ + 2se''^ + se-''^ + s^e''^ 

3 3 P + SPs'^e'^''^ + 3se''^ + s^e'^''^ 

4 2 cos {2ip) + 2 (2 Ps^ + 2 Ps) cos (^p) + 1 + 4 + 

4 2 l + 3Ps3e'*' + Ps^e^''^+3Pse'^ + Pse-'*' + 3s2+3s^e^''^ + s*e^''^ 

4 4 l + 4Ps^e3^'^+4Pse^'^ + 6s2e"^ + s*e'''^ 

5 1 P + 2 Ps*e^ + 6 Ps^ + 6 s^e^*" + 3 s^e"^"^ + 2 se"'"^ + s^e'*^ + Ps^e'^ + 3 Ps'' + 3 Ps^e^ '"^ + 3 se*"^ + s^e^ '"^ 
5 3 P + 4 Ps^e^ + 4 Ps^ + 6 s^e'"^ + se-''^ + s^e^ + Ps'^e* '"^ + Q Ps^e^ + 4 se'-^ + 4 s^e^ '"^ 

5 5 P + s^e^ + 5 Ps^e* '-^ + 10 Ps2e2'^ + 5 se'^ + 10 8^6^'"^ 

6 2 Ps^ cos {3(p) + 2 (3 + 3 s*) cos (2 v?) + 2 (3 Ps + 3 Ps^ + 9 Ps'^) cos (v?) + 1 + 9 s* + 9 + s*^ 

6 2 1 + 12 Ps^e^-^ + 4 Ps^e"'^ + 4 Ps^e^ + 4 Pse'"^ + 2 Pse"'^ + 8 + 2 Ps^e^ + s'^e^ + 4 Ps^'^e''^ + s^e'^ '"^ + 

6 s^e"^ + 8 s*e2 + s*e* + 6 

6 4 1 + W Ps'^e^'^ + 10 Ps^e^^'^ + 5 Pse''^ + Pse~''^ + 5 +5 Ps^e^^'^ + Ps^'e^''^ + s^e*^'^ + 10 s^e'^''^ + 10 s'^e'^''^ + 

6 6 l + 20Ps^e^'^+6Pse'^ + 6Ps5e5^'^ + s''e''^'^ + 15s2e2^'^ + 15/e*''^ 

7 1 P + 12 Ps^e^ ''^ + 12 Ps^ + 4 Ps^e'^ '"^ + 18 s^e'^ + 12 s^e"'^ + s^e"^ '"^ + 3 se"'"^ + 3 s^e^ '"^ + 12 s'^e^'^ + 6 s^e"^"^ + 

Ps^e*^"^ + 3 Ps^e^^ '"^ + 18 Ps^ + 6 Ps^e^ '"^ + 4 se"*' + 4 s^e^ + 3 Ps^'e^ + s'^e''^ + 4 Ps^ 

7 3 P + 20 Ps''e2'^ + 10 Ps^ + 2 Ps^'e''^'^ + 20 8^6'"^ + 5 s^'e-''^ + 2 se~"'^ + 10 3^6^^'^ + 10 s^e^-^ + 5 Ps*e*''^ + Ps2e-2^^ 

10 P/ + 10 Ps^e^ + s'e^ '"^ + 5 se'"^ + s^^e^ + 10 s^e^ + 5 Ps^e^ '"^ 

7 5 P + 20 Ps^e^ + 6 Ps^ + 6 Ps^'e'' + Ps^e^ ''^ + 15 s^e^*' + se""*' + 15 s^e^ ''^ + 15 Ps^e'' + 15 Ps^e^ + s' + 

7 7 P + 7Ps''e«"^ +35 P/e''''^ + 21 Ps^e^'*' + sV"^ +7se"^ + 21 s^e'^*'^ +35 8^6^''^ 

TABLE II: The primitive basis function in the integral (I47p translated to {a, ip) coordinates. The cases of negative m follow 
immediately by complex-conjugation. 



where tan t^o- = <yyl <^x defines the azimuth of the wave number. The Bessel Function integrals are recursively computed 
via [il, 11.3.4] 

= t J,n{t)dt ^ ^ y (52) 

Jo \ a"^'^J„i+i[a) + [m ~ n)g(n,m + l,a) ,n>m. 

The integrals j4„_„i(g) in ^ are just the special case of zero momentum, ct = 0. In this limit, (|50|) simplifies to 
Is<q exp{i'irup)(Ps ~ 27r(5,„,o(7"^^/ + 2), which proposes an alternative to compute Table [H 

We visualize the Fourier representations of the first 12 basis functions in Figures [5] and [Sj The fj{x,y) have a well 
defined parity with respect to refiection across the origin, i.e., the terms in (p7|) - ((39)) are superpositions of z" sin(m6') 
and z" cos(m0) with a fixed, common parity (—1)™ in each line. As a consequence, the fj{crx,cry) are purely real or 
purely imaginary, and only the non- vanishing of the two components is shown. 

/i in Figure [S] is the familiar diffraction limited point spread function of the two-beam interferometer [l6j . /i in 
Figure [2] is the convolution of the double pinhole mask with a single circular telescope pupil of area irq^ . Its Fourier 
representation is the product of the circular Airy disk centered at cr = by the hyperbolic pattern of fringes with 
their narrower width along as determined by the baseline length. 



VII. SUMMARY 



Orthogonality of functions over the area of two non-overlapping circles is defined according to the algebraic stan- 
dards. The fundamental areal integral for functions that are of simple analytical format in the spherical coordinate 
system attached to the center of symmetry has been recursively reduced to polynomials or Gaussian hypergeometric 
functions of the normalized circle radius. This defines a set of orthogonal basis functions over the common area, if 
the Zernike basis functions, for example, provide the input to the Gram-Schmidt procedure. The first few of these 
have been written down in analytical form. For the higher-order aberrations, the complete information for numerical 
instantiation has been presented. 
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FIG. 5: The Fourier transforms of basis functions /i to fe for q = 1/2. 



On a side note, computation of the interferometric signal of the pupil-beam recombination given an expansion in 
the radial-azimuthal coordinates has been reduced to a double sum of over a generic overlap integral. 



APPENDIX A: AUXILIARY INTEGRAL B 



1. Two Recursion Strategies. 



This section deals with the evaluation of the integrals defined in Splitting off and expanding a square of the 

(n + 2)nd power, we obtain 



where we have introduced 



arcsin q 



cost 



cos 9 ± \/ q"^ — sin 



n+l 



cos{m9)d9. 



(Al) 



(A2) 



This defines a first strategy to evaluate Bn+2,m{q)' recursive reduction of the first lower index in steps of two at the 
expense of implementing the Xn+i^mil) integrals — those to be treated in Appendix IB] Starting the recursion from 
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FIG. 6: The Fourier transforms of basis functions fr to /12 for q — 1/2. 



an even n, one ends up at 



0,m 



(q) = — sin(m arcsing) 
m 



(A3) 



This is also correct in the hmit m — > 0, and equivalent to zero as the difference Bq^^ — Bq^^ is formed. The symbol ~ 
is reserved in this script to indicate that terms on the right hand side have been removed which cancel if differences 
B~^{q) — B~{q) (common pair of subscripts with both B) or differences X^{q) — X^{q) (common pair of subscripts 
with both X) are calculated. For odd n, the recursion (|Aip terminates at 



BtrM 



dO cos 6 ± y q"^ — sin 6 ) cos{m9) 



± 



arcsm q 



dOJq'^ - sin-" e'T„(cos( 



(A4) 
(A5) 



For m — this value will be given in (jA16p . For larger m, a mixture of partial integrations and the product rule for 
Chebyshev polynomials [J, 22.7.24] generates a recursion for the second index, again in steps of two |18l |: 



R± . . , 2(m-2)(l-2g^) ^ 



2(9) 



m + 1 



(A6) 
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Values to start this recursion at small even m are discussed in Section lA 21 — this implies that n + m is odd and does 
not happen for the Gram-Schmidt seeds proposed in Section HVl At small odd m the recursion starts from 



B^.iq) = ±-g2; B^.iq) = ±-q'{l - q'). (A7) 

An alternative second strategy to evaluate Bn-\.2,mi<l) looks as follows: Binomial expansion of the (n + 2)nd power 
in the integrand of (fT4|) yields 

L(n+1)/2J «+l-2s . , „ 

'n — s + l\, o v-^' / n + 1 — 2s 



Bt^2^m{<l) - ± E H'^' - 1)^ E U±li_ Bt+rrM + Bt,^_^,iq) • (A8) 

n—u odd 

Here, the prime at the sum symbol means the term for u = 0, if it occurs, is to be halved. The differences to the first 
strategy are 

• No evaluation of the X^^^ mil) needed. The contents of Appendix [Bl can be ignored. 

• The decrement of the first index of B^^2 m (|A8[) down to 1 comes at the cost of an increment of some second 



indices, eventually a recourse to (|A6[) . 

Use of (jA8[) in conjunction with (jA6[) eventually uses (|A7[) . whereas (|A7[) and the case of odd m are made 
irrelevant for the first strategy through (fT5)l . 



Explicit examples of (jASp for small n are: 

Bf^^ ~ Bi^m+i + m > 2; (A9) 

Bt,m — + Bi,„_2 + (l + g')Si,™, m>2; (AlO) 

5±2 ~ 2(1 + q^)Bi,^ + (1 + 2(72)^1,3 + -Bi,5; (All) 

Bf .^ ~ Bi,„_3 + (1 + 2g2)Bi,„_i + (1 + 2g2)Bi,„+i + Bi,„+3, m > 3; (A12) 

S5±2 ~ (1 + 3g2)Bi,o + (2 + V + 5^ 2 + (1 + 39')Bi,4 + Bi,6; (A13) 

B5±3 ~ (2 + 3g2) + (1 + Aq^ + q^) Si.3 + (1 + 3g2)Bi.5 + Bi,7; (A14) 

B5±„ ~ Bi,™_4 + (l + 3g2)Bi,„,_2+(l+V+g4)5^^^ (A15) 

+ (1 + 3g2)Si,,„+2 + Si,™+4, m > 4. 
The argument q at all the B, ,{q) has been omitted for brevity. 

2. The case of odd n + m. 

With a substitution sin0 — ^, i?i^o can be written as a superposition of complete Elliptic Integrals of the first and 
second kind: 



ccos yjl — cp 

Bfjq) = / (cos0± Vg2-sin^ 







q ± / ^^^ rfC = q ± [£^(9^) - (1 - g')i^(«')]- (A16) 
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Bt^q) ^ ME{q') - (1 - q')K{q')]. (A17) 
A merger of the series expansions of the EUiptic Integrals yields 

qy m-m? _ 1/2,1/2 



^1,0(9) ^ I'^V2J -n(7TT)r" 4'^ \ 2 



25 245 o 1323 ,„ 

1024^ 16384^ 131072^ 



7623 , 184041 _^ 4601025 , 29548805 , 193947611 ^20 



9-1 



1048576" 33554432" 1073741824" 8589934592" 68719476736 
±1. 



Due to a logarithmic singularity at — 1, the power series converges poorly if the argument of the Elliptic Integrals 
approaches unity [l^l ■ 

(jA17p is the first anchor value for (|A6[) . The second is 



/■arcsin q y 

B^^iq) - ± / Jq^-sm^qcos{2e)de 
Jo 



(A18) 



±1(1 - ± 1(2^2 _ 1)^(^2)^ ^A19) 

With the aid of (|A6p . expressions for B^^-^(q) are bootstrapped from (|A17p and (|A19p . This sequence starts: 

BfM ^ ±-5-(g2 _ 1)(8q2 _ ± i-(-16g4 + iQq^ - l)E{q^)- (A20) 

15 15 

BtM) - - l)(128g4 - 80g2 + i)K{q^) ± ^(2g2 _ l)(128g4 - I28g2 + i)E{q^)- (A21) 

Bfjq) ~ ± — (g2-i)(i0249^-1152/ + 288g2-5)/s:(g^)± — (-2048g** + 4096g^-2496/ + 448q^-5)i;((7^). (A22) 
315 315 



APPENDIX B: AUXILIARY INTEGRAL X 



The auxiliary integrals (|A2p are put into an algebraic format by the substitution sin0 = qz, 

xt+i,m{q) = q {V^^l^±qVT^y'''TUVi^l^)dz, (Bi) 

and then broken down through binomial expansion of the (n + l)st power and the explicit polynomial expression for 
the Chebyshev function 1, 22.3.6] via 9, 3.197.3] 



L«/2J 



X^^,^^{q)^±r--\{n + iy.q' ^ 



1 



3=0 



in-2jy.j\ \2 



2j 



Lm/2J/ Iw 



J2 ( 



'4) 2(m-cr) 



X < 



^rr)^r ;™>o; 



F I ^ 2' 2 I q2 



: TO = 0. 



If n + TO is an even number — which is the case for the Gram-Schmidt procedure described in Section IIV BI — the 
hypergeometric series terminate and become polynomials of q^ of order (n + to)/2 — j — a [ll, 15.4.1]. The result can 
be tabulated in terms of power series coefficients aj{n, to), 



^n+i,miq) = T^q^Oij{n,m){±qy. 

J=0 



(B2) 



Only the values with odd j are of interest, because our application eventually looks only at the differences X^^-^ miq) ' 
■^n+i miq) i'^ which terms of even j cancel. The basic values for m = are in Table Hill 



17 



i {n, 0) 

1 (ri + l)/4 

3 n(n + l)(n - 2)/32 

5 n(n + l)(n-2)2(n-4)/768 

7 n(n + l)(n-2)2(n-4)^(n- 6)/36864 

9 n(n + l)(n - 2)2(n - 4f{n - 6)^(n - 8)/2949120 



11 n(n + l)(n - 2)2(n - 4)^(n - 6)^(n - 8)^(n - 10)/353894400 
TABLE III; Expansion coefficients of (|B2]| at m = 0. 



Table [TVl summarizes (IB2[) for small values of n and m. 

n m T,j^i,3,5,7,...<^An,m)q^ 

+1/4(7 

2 +1/4(7 - l/8g^ 

4 +l/4g - 1/2(7^ + l/4<j5 

6 +1/4(7 -9/8g3 +3/29^5 -5/8g'^ 

8 +l/-iq-2q^ + 5q^ -5q'' + 7/iq^ 

1 +l/2g- 1/16(7^ - 1/128(7-^ - 5/2048(7'^ - 35/32768g^ - 147/262144g" 

-693/2097152^^=* - 14157/67108864^^^ - 306735/2147483648g" + . . . 
1 2 +1/2(7 -5/16(7^ + 7/128(7^ + 15/2048g'^+ 77/32768,7^ +273/262144,7" 

+1155/2097152(7^=* + 21879/67108864?^^ +448305/2147483648?^^ + . . . 
1 4 +1/2(7 - 17/16g=* + 95/128g^ - 245/2048,7^ - 483/32768(7^ 

-1155/262144(7" - 3861/2097152g^=* - 62205/67108864(7^^ 

-1130415/2147483648?" + . . . 
1 6 +l/2g - 37/16(7=* + 455/128?^ - 4305/2048?'' + 9933/32768?^ 

+9009/262144?" + 20163/2097152?^=* + 255255/67108864?^^ 

+3926065/2147483648?" + . . . 

1 8 +1/2? - 65/16?=* + 1407/128?^ - 26565/2048?^ + 213213/32768?^ 

-219219/262144?" - 182325/2097152?" - 1524237/67108864?" 
-18244655/2147483648?" + . . . 

2 +3/4? 

2 2 +3/4? -3/8?=* + 1/8?=^ 

2 4 +3/4?-3/2?=* + 5/4?5-3/8?'^ 

2 6 +3/4? - 27/8?=* + 45/8?^ - 33/8?^^ + 9/8?^ 

2 8 +3/4? - 6?=* + 17?-'^ - 45/2?^^ + 57/4?^ - 7/2?" 

3 +? + 3/8?=* - 1/64?^ - 1/1024?^^ - 3/16384?^ - 7/131072?" 

-21/1048576?" - 297/33554432?" - 4719/1073741824?" + . . . 
3 2 +? - 1/8?=* + 7/64?-""' - 13/1024?^^ - 19/16384?^* - 35/131072?" 

-93/1048576?" - 1221/33554432?" - 18447/1073741824?" + . . . 
3 4 +? - 13/8?=* + 95/64?=^ - 625/1024?^ + 1085/16384?^ 

+777/131072?" + 1419/1048576?" + 15015/33554432?" 

+196625/1073741824?" + . . . 
3 6 +? - 33/8?=* + 455/64?=^ - 6125/1024?^ + 37485/16384?^ 

-30723/131072?" - 21021/1048576?" - 148005/33554432?" 

-1519375/1073741824?" + . . . 

3 8 +? - 61/8?=* + 1407/64?=^ - 32193/1024?^ + 383229/16384?'-' 

-1072071/131072?" +819819/1048576?" +2122263/33554432?" 
+ 14272401/1073741824?" + . . . 

4 +5/4? + 5/4?=* 
4 2 +5/4? + 5/8?=* 

4 4 +5/4?-5/4?=* + 5/4?=5-5/8?'' + l/8?^ 

4 6 +5/4? - 35/8?=* + 15/2?^ - 55/8?^" + 13/4?^ - 5/8?" 

4 8 +5/4? - 35/4?=* + 25?'^ - 75/2?^ + 125/4?'' -55/4?" +5/2?" 

5 +3/2? + 45/16?=* + 45/128?^ - 15/2048?'^ - 9/32768?^ - 9/262144?" 

-15/2097152?" - 135/67108864?" - 1485/2147483648?" + . . . 
5 2 +3/2? + 33/16?=* +5/128?^ +45/2048?^ - 57/32768?^ - 29/262144?" 
-39/2097152?" - 315/67108864?" - 3245/2147483648?" + . . . 



+3/2g - 3/16(j^ + 77/128g^ - 735/2048g'^ + 3255/32768?'' 
-2009/262144(7" - 1071/2097152?^^ - 6039/67108864?^'^ 
-50765/2147483648?" + . . . 

+3/2? - 63/16?'^ + 837/128?^^ - 12915/2048?'' + 112455/32768?'' 
-242109/262144?" + 147609/2097152?" +312741/67108864?^^ 

+1756755/2147483648?" + . . . 

+3/2? - 147/16?=^ + 3245/128?^^ - 79695/2048?^" + 1149687/32768?" 
-4749129/262144?" + 9810801/2097152?" - 23146695/67108864?" 
-47732685/2147483648?" + . . . 
+7/4? + 21/4?3 + 7/4?5 
+7/4? + 35/8?=^ + 7/8?^ 
+7/4? + 7/4?3 

+7/4? - 21/8?^ + 35/8?-' - 35/8?"^ + 21/8?" - 7/8?" + l/8?^''* 
+7/4? - 35/4?=^ + 91/4?^ - 35?^ + 133/4?" - 77/4?" + 25/4?" 
-7/8?" 

+2? + 35/4?^^ + 175/32?^^ + 175/512?"^ - 35/8192?" - 7/65536?" 

-5/524288?^-^ - 25/16777216?^-''' - 175/536870912?" + . . . 

+2? + 31/4?'^ + 119/32?""' + 35/512?"^ + 77/8192?" - 35/65536?" 

-13/524288?" - 53/16777216?" - 335/536870912?" + . . . 

+2? + 19/4?^ + 15/32?^ + 63/512?'" - 483/8192?" + 777/65536?" 

-357/524288?" - 585/16777216?" - 2607/536870912?" + . . . 

+2? - 1/4?^^ + 55/32?""' - 957/512?^ + 9933/8192?" - 30723/65536?" 

+49203/524288?" - 92565/16777216?^" - 155727/536870912?" + . . . 

+2? - 29/4?^^ + 559/32?^^ - 13585/512?^ + 213213/8192?" 

-1072071/65536?" +3270267/524288?" -20796633/16777216?" 

+38955345/536870912?" + . . . 

+9/4? + 27/2?^ + 27/2?^ + 9/4?^ 

+9/4? + 99/8?-^ + 21/2?^ + 9/8?^" 

+9/4? + 9?^ + 15/4?5 

+9/4? + 27/8?^^ 

+9/4? - 9/2?=* + 21/2?^ - 63/4?^" + 63/4?" - 21/2?" + 9/2?" 

-9/8?" + 1/8?" 

TABLE IV: Table of aj{n,m)q^ , summed over odd j only. For odd 
n + m, the series is shown up to 0{q^^), indicated by the triple dots. 
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APPENDIX C: NOTATIONS 

[J floor function; largest integer not greater than the argument 

(., .) scalar (inner) product between the two arguments pB]) 

"Yl! u= summation with the term of index m = halved 

~ equivalent upon subtraction of or of superscripts of opposite sign 

An,m{q) the integral (O 

ctn,m{<l) the expansion of the X integral (|B2p 

B^.^{q) the integral dH]) 

B{x, y) = r(a;)r(y)/r(a; + y), the Beta Integral [9, 8.38] 

5jk Kronecker delta; equal to \ li j = fc, equal to if 7^ fc 

E{.) complete Elhptic Integral of the first kind [1, §17.3] 

E electric field amplitude 

fk{q) fcth orthogonal basis function of the binocular area 

F( '^'^ I .) Gaussian hypergeometric function [1, §15] 
c 

azimuth angle in spherical coordinates centered at circular sub-pupil 

7._. abscissa sections of the Gram-Schmidt procedure 

gk fcth input ("seed") function to the Gram-Schmidt synthesis 

(7.^. Bessel Function integral [1_, 11.3.1] 

/ interferometric intensity 

i imaginary unit 

J(.) Bessel function of the First Kind 

j\\ double factorial, = 1 • 3 • 5 • 7 • • • j if j is odd, = 2 • 4 • 6 • • • j if j even. (—1)!! = 1. 

K{.) complete Elliptic Integral of the second kind ]1, §17.3] 

q circle diameter in units of circle center distance 

s radial distance to a circle center, < s < q 

(T, (jx, (Jy, a wave number, Cartesian coordinates, modulus 

Tm{-) Chebyshev Polynomial the first kind of order m [i, §22] 

azimuth angle in the circular coordinates centered at midpoint between sub-pupils 

^trrM) the integral (lA2|) 

Zk{.) Zcrnike circle functions [17] 

z radial distance to a center of symmetry < z < 1 + g 
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